Non-tumor-related prognostic factors for immunotherapy–chemotherapy or immunotherapy alone as first-line in advanced non-small cell lung cancer (NSCLC)

Besides programmed death ligand 1 (PD-L1) expression, rapid, cost-effective and validated scores or models are critical for the prognosis and prediction of patients received immune checkpoint inhibitors (ICIs). In this retrospective study, 182 patients with NSCLC receiving ICIs from 2015 to 2022 were divided 1:1 into a training cohort and a validation cohort. We identified a score established by three factors and analyzed the prognostic implications by Kaplan–Meier approach (Log rank test) and time-dependent receiver operating characteristic (ROC) analyses. A non-tumor-related score (NTRS) was established that could be used as a prognostic factor (HR 2.260, 95% CI 1.559–3.276, P < 0.001 in training cohort; HR 2.114, 95% CI 1.493–2.994, P < 0.001 in validation cohort) and had a high time-dependent ROC for overall survival (OS) (AUC 0.670–0.782 in training cohort; AUC 0.682–0.841 in validation cohort). PD-L1 (1–49%) and NTRS (score = 0, 1, 2, 3) combination significantly improved the assessment of patients’ OS and progress-free survival (PFS), which was statistically different in training cohorts (P < 0.001 for OS, 0.012 for PFS) and validation cohorts (P = 0.01 for OS, < 0.001 for PFS). The NTRS provided a better assessment of durable clinical benefit (DCB) compared to PD-L1 expression (P = 0.009 vs. 0.232 in training cohort; P = 0.004 vs. 0.434 in validation cohort). NTRS may help improve prognosis stratification of patients receiving ICIs in first-line NSCLC and may be combined with tumor-related parameters. Supplementary Information The online version contains supplementary material available at 10.1007/s10238-024-01298-z.


Introduction
ICIs have revolutionized the treatment landscape of NSCLC.ICIs have been approved for first-line therapy and are available either alone (PD-L1 tumor percentage score (TPS) ≧ 50%) or in combination with chemotherapy [1,2].Several prognostic and predictive factors in NSCLC have been described in association with ICIs.Some clinical studies have shown that PD-L1, tumor mutational burden (TMB), microsatellite instability (MSI), defective mismatch repair (dMMR) and other immune features are considered as potential biomarkers for predicting ICIs [3][4][5][6].However, only high levels of MSI/dMMR and TMB (no universal definition of high TMB; high cost; limited by tissue sample size) are expressed in a minority of patients [3,[7][8][9].Despite Keynote-024 [10] and Keynote-042 [11] demonstrated a benefit in OS in patients receiving ICIs with PD-L1 ≧ 50%, patients in Keynote-189 [12] and 407 [13] with ICIs combined with chemotherapy, regardless of PD-L1, had a significantly superior objective response rate (ORR), PFS, and OS.Therefore, the predictive and prognostic value of PD-L1 expression is limited.
Despite the tremendous progress of ICIs in the treatment of NSCLC, there are still many outstanding issues in patient prognosis and efficacy that need to be solved.Some studies showed genetic models [14,15], neural network model [16], and nutrition index [17] could predict the risk of NSCLC progression.However, many models had not reached clinical practice.Obviously, this was related to economic and methodological complexity and/or tumor heterogeneity.There is growing knowledge about the importance of comprehensive prediction as a biomarker or model.Clinically relevant prognostic models are also being performed in full speed.Appropriate clinical models in early NSCLC with adjuvant chemotherapy [18], advanced NSCLC in old age [19] and lymph node metastases (LNMs) [17] have been studied as prognostic models.Nevertheless, non-tumor-related predictive and prognostic models without specific population restrictions in advanced NSCLC patients receiving ICIs are controversial.Worse yet, the relative importance of prognostic factors in the era of ICIs is currently not well defined and explored.In the clinical field, there is an urgent need for valid, simple, cost-effective, and comprehensive model predictions to identify patients who may benefit from ICIs.
This study aimed to explore potential model and to accurately predict prognosis in patients with advanced NSCLC.A new NTRS could more precisely assess the predictive and prognostic value of NSCLC patients undergoing ICIs and more clearly identify the relative prognostic impact of clinical factors.The NTRS established a new approach that can be applied to the prediction and prognosis of NSCLC.

Materials and methods
This retrospective study was approved by the Research Ethic Committee of Peking University International Hospital's, and all patients provided their informed consent to treatment including immunotherapy.The process followed the declaration of Helsinki.

Study population
We retrospectively screened 1796 NSCLC patients at XX Hospital from Jan. 2015 to Jun. 2022.All patients diagnosed with histologically confirmed stage IV or unresectable and non-curative stage III NSCLC who received at least two doses of first-line therapy were enrolled.Inclusion criteria were histologically or cytologically confirmed NSCLC with at least one measurable tumor lesion according to response evaluation criteria in solid tumors (RECIST) version 1.1, adequate bone marrow and organ function, and an expected survival time of ≧ 3 months.Exclusion criteria were patients with tumors with targetable mutations (e.g., EGFR mutations, ROS1 and ALK translocations), patients enrolled in clinical trials, patients receiving antitumor therapy such as chemotherapy, radiation, targeted and immunotherapy prior to the first dose, use of systemic glucocorticoids 4 weeks prior to the first dose, and metabolic syndrome such as type 2 diabetes, hypertension, family history of obesity, atherosclerosis, fatty liver and BMI > 27.Due to patient blood tests resulted after diagnosis, patients were placed on a relatively standardized diet with a reference daily sugar, protein and lipid ratio of 3:2:1.Enrollment data included patient characteristics (age, gender, smoking history, Eastern Cooperative Oncology Group (ECOG) status, blood tests, biochemical tests, coagulation tests and tumor marker tests), tumor characteristics (histology, stage, metastatic sites, and PD-L1 expression, baseline tumor size (BTS), best overall response (BOR), and DCB).

Blood sample analysis
Blood tests performed within 1 month prior to treatment were involved in the analysis, and laboratory blood tests were performed.The SYSMEX XN3000 analyzer was used to determine the levels of routine blood (white blood cells, lymphocytes, neutrophils, platelets and hematocrit (HCT)).The BECKMAN AU5800 automatic biochemical analyzer measured lactate dehydrogenase (LDH) (the upper limit of normal (ULN) of 250 U/L), alkaline phosphatase (ALP) (ULN of 125 U/L), albumin (ALB) (the lower limit of normal (LLN) of 40 g/L), high-density lipoprotein cholesterol (HDL-C) (LLN of 1.0 mmol/L), low-density lipoprotein cholesterol (LDL-C) (ULN of 3.4 mmol/L), free fatty acids (ULN of 769 μEq/L), and β2-microglobulin (ULN of 2.64 mg/L).ACL-TOP 700 automatic coagulation analyzer was used to test plasma fibrinogen (FIB) (ULN of 400 ng/ dL).CobasE601 analyzer examined carcinoembryonic antigen (CEA) (ULN of 5 ng/mL), squamous epithelial cell carcinoma-associated antigens (SCC) (ULN of 1.5 ng/mL), and CYFRA21-1 (ULN of 3.3 ng/mL).Each patient underwent baseline chest and abdominopelvic computed tomography (CT) scans as well as magnetic resonance imaging (MRI) of the brain and bone scans as necessary less than 8 weeks prior to the first ICI.

Scoring model development and validation
They were randomized splitting into training and validation groups and then stratified by gender and age.NTRS was developed in the training cohort.To determine the prognostic significance of each factor in the baseline clinical variables (age, gender, smoking history, pathology type, ECOG status, CEA, SCC, CYFRA21-1, size, T-stage, N-stage, metastatic sites, and PD-L1 expression), the best cutoff values were screened and validated to confirm the prognostic ability.Three significant factors were eventually included in the score.One point was scored for each of these three variables, and the total score was defined as the scoring model.The discriminative ability of the model was assessed by the area under the curve (AUC).The accuracy of the model was determined by the Harrell's Discriminant Consistency Index (C-index).

Efficacy assessment and treatment modalities
These measurable lesions were performed according to RECIST.Based on National Comprehensive Cancer Network (NCCN) guideline recommendations and PD-L1 expression, patients were treated with single-agent immunization or immune-chemotherapy treatment.In addition, patients with ECOG≧ 2 partially received immune-monotherapy according to family wishes.

Follow-up
Patients were followed up primarily by outpatient visits (or telephone) every 10 ± 2 weeks.PFS was defined as the time from initiation of treatment to disease progression or death.OS was defined as the time from initiation of treatment to death.DCB was determined as a complete response (CR), partial response (PR) or stable disease (SD) lasting longer than 6 months; no durable benefit (NDB) was characterized as progressed disease (PD) or SD lasting less than 6 months.BOR was measured as CR/SD/PR groups.Patients with no events were censored at the last follow-up.Missing data rates were < 15% to 20%.

Statistical analysis
IBM SPSS (version 26.0; SPSS Inc, Chicago, IL, USA), MedCalc (version 20.1.0)and R software (http:// www.rproje ct.org/) version 3.6.1 were used to conduct statistical analysis.Continuous variables were presented by median with interquartile ranges (IQRs) and categorical variables were expressed as frequencies (%).Continuous and categorical variables were assessed by two-side Mann-Whitney U test and chi-square test or Fisher's exact test, respectively.Kaplan-Meier method was used to estimate PFS and OS, and differences were compared using the log-rank test.X-tile [23] analysis was determined as an identification of the optimal cutoff for 5-year OS with factors including neutrophil, lymphocyte, high neutrophil lymphocyte ratio (NLR), platelets, platelets lymphocyte ratio (PLR), HCT, β2-microglobulin, FIB, LDH, ALB, HDL-C, ALP, CYFRA21-1, SCC, CEA.Univariable and multivariable Cox regression analyses were performed to determine independent factors of OS or PFS by the Cox regression methods.Time-dependent ROC curve analyses used the "timeROC" package to evaluate the prognostic and predictive ability of each factor.P values < 0.05 were considered statistically different.

Identification of the NTRS with OS in the training and validation cohort
Univariate and multivariate Cox regression analyses in Table 2 and Table S1 demonstrated the relationship between clinicopathological variables and OS and PFS in the training and validation cohorts.In the univariate analysis of the training and validation cohorts, PD-L1 expression, maximum diameter, LDH levels, CYFRA21-1 levels and NTRS were associated with OS; NTRS was associated with PFS.These variables were subsequently included in multivariate analyses.
Thus, even though PD-L1 had statistically different OS in the < 1%, 1-49%, and ≧ 50% groups, there was no PFS difference.Moreover, both OS and PFS were statistically different in the PD-L1 1-49% stratification in each NRTS score (0-3) in both the training and validation cohorts.

NTRS showed higher efficiency than PD-L1 in predicting response
NTRS demonstrated a higher efficiency than PD-L1 expression as a biomarker for immunotherapy.When response was evaluated with RECIST1.1,NTRS were more different between the DCB and NDB groups than PD-L1 expression (Fig. 4A-D), but there was no significant difference between the BOR and non-BOR groups (Fig. 4E-H).The proportion of NTRS and PD-L1 was 90% (N = 81) and 90% (N = 81) for non-PD (PR + SD) patients in the training cohort and 92% (N = 83) and 91% (N = 82) for non-PD patients in the validation cohort (Fig. 4I-L).
ROC analysis demonstrated that NTRS significantly enhanced the predictive efficiency of PD-L1 expression (Fig. S5).When BOR was assessed in the training cohort, the AUC of NTRS was 0.650, which was the trend of assessing prognosis than 0.522 of PD-L1 (Figs.S5A, 5E).In evaluation of DCB in the training cohort, the AUC of NTRS was higher, 0.644, than 0.572 for PD-L1 (Fig. S5B,  5E).When BOR was assessed in the validation cohort, the AUC for NTRS was 0.584, significantly lower than 0.637 for PD-L1 (Figs.S5C, 5E).In evaluation of DCB in the validation cohort, the AUC of NTRS was highest at 0.665, which was higher than that of 0.547 for PD-L1 (Fig. S5D, 5E).CI, confidence interval; OS, overall survival; PFS, progression-free survival; IQR, interquartile range; ECOG, Eastern Cooperative Oncology Group; CEA, carcinoembryonic antigen; SCC, squamous epithelial cell carcinoma-associated antigens; LDH, lactate dehydrogenase; ULN, upper limit of normal; ALP, alkaline phosphatase; LUAD, lung adenocarcinoma; LUSC, lung squamous cell carcinoma; BTS, baseline tumor size; PD-L1, programmed cell death ligand 1; TPS, tumor cell proportion score Taken together, the combination of PD-L1 expression and NTRS had a greater AUC and higher predictive prognosis value to improve the effectiveness of PD-L1 expression as an immunotherapy biomarker.

Distribution and classification characteristics of NTRS and PD-L1 expression in metastatic sites and tumorigenic factors
The distribution of PD-L1 expression in tissue samples from different metastatic sites was moderately different depending on organ differences (Fig. S6).Among metastatic sites, the proportion of patients with high PD-L1 expression was greatest in lung (PD-L1 (≧ 50%) = 29%) and smallest in pleural (PD-L1 (≧ 50%) = 13%), and the proportion of patients was greatest in the adrenals (scores (1 + 2 + 3) of 66%) and smallest in the lung (scores (1 + 2 + 3) of 29%) (Fig. S6A, B).Differences in the distribution of PD-L1 expression and NTRS at different metastatic sites were further explored.PD-L1 status (PD-L1≧ 1% vs. < 1%) was predictive of ICI treatment response in brain metastases, but no significance was found between metastasis site and treatment response in the risk score (Fig. S6C, D).Categorical and regression analyses were used to assess the prognostic factors of NTRS and PD-L1 expression in long-term survival.Table S2 provided the tumor correlates of NTRS and PD-L1 expression.In PD-L1 expression, the level of the tumor marker SCC was identified as the most important correlate, while in the NTRS, LNMs status was determined as the most important correlate.And the critical ranges of LNMs and SCC levels were identified as N2 and 2.5 ng/ml.Specifically, regression analysis of critical values in the PD-L1≧ 50% group suggested worse PFS when SCC levels exceeded 2.5 ng/ml (mPFS: 8.4 vs. 2.4 months, P < 0.001), but no difference was found in OS (Fig. S7).In addition, in the NTRS, LNMs suggested survival differences in OS and PFS (Fig. 5).In particular, N2 staging resulted in shorter OS and PFS in this group of patients with NTRS ≧ 1 (mOS: 42.6 vs. 14.2 months, P < 0.001; mPFS: 10.37 vs. 4.17 months, P = 0.001).These results implied longer OS in patients with PD-L1 ≥ 1% brain metastases, but there was no correlation between NRTS and metastatic site.In a further, patients with NTRS (≧ 1, host factor) and LNMs (N2, tumorigenic factor) had worse OS and PFS; patients with PD-L1 (≧ 50%) and SCC levels (> 2.5 ng/mL) had shorter PFS.

Subgroup analyses of OS and PFS based on NTRS and PD-L1 expression
Patients with NTRS scores (0 vs. ≧ 1) (Fig. 6) and PD-L1 expression (0 vs. ≧ 1) (Fig. S8) were examined for differences in OS and PFS in the training and validation cohorts, respectively.When risk scores ≧ 1, OS was significantly shorter for both risk factors (female: P = 0.01; size < 86 mm: P < 0.025) (Fig. 6A, C), but PFS was not uniform for these factors in the training and validation cohorts (Fig. 6B, D).Furthermore, in the five factors of male, size < 86 mm, ECOG 0-1, N2-3, and SCC level less than 2.5, OS was statistically prolonged despite the absence of positive expression of PD-L1 in the training and validation cohorts (Fig. S8).To summarize, OS shortened with increasing PD-L1 levels and NTRS, but no similar results for PFS.Also, there are 5 factors (male, size < 86 mm, ECOG 0-1, N2-3, SCC level < 2.5) that indicated an extended OS.

Discussion
ICIs have been fundamentally altering the administration of NSCLC.Despite the success of ICIs in a proportion of NSCLC patients, some patients continue to progress and die in a relatively short-term period.A major challenge is to establish which patients will respond beneficially to ICIs and have a longer prognosis.This paper has provided an alternative foundation for the analysis of the prognostic and predictive value of NTRS that could improve more significant survival and benefit from ICIs in a fuzzy environment.
Emerging studies of risk models [14][15][16][17][18][19] have been introduced to predict the risk of NSCLC progression.The survival neural network model investigated more promising results in the prediction of lung cancer-specific survival than the tumor, node, and metastasis stage (C statistic = 0.739 vs. 0.706) [16].According to a model based on ferroptosisrelated genes, a series of bioinformatics analyses were conducted to verify the predictive ability of OS (P < 0.001 vs. 0.0095, AUC = 0.679 vs. 0.618 in ≧ high and low groups between training and validation cohort) as independent prognostic factors (HR 0.63, 95% CI 0.44-0.91,P = 0.01) for NSCLC patients with immunosuppressive status [14].High tumor microenvironment related genes model was identified to predict longer OS (log-rank P = 0.044) and was an independent factor (HR 3.32, 95% CI 2.16-5.09)[15].CT-based imaging signs, systemic inflammatory response index, prognostic nutritional index and CEA were filtered and included in the nomogram (AUC 0.828, 95% CI 0.778-0.877 in the training cohort, AUC 0.816, 95% CI 0.737-0.895 in the validation cohort, respectively) for predicting LNMs in NSCLC patients [17].Patients were classified into good (0-8) and poor (≧ 9) groups based on clinical prognostic model (3 (if male) + 3 (if ECOG = 1) + 8 (if ECOG = 2) + 11 (if initial stage = IV) + 4 (if weight loss)) [19].The median OS in the two groups were 13.15 (95% CI 10.82-15.91)and 8.52 months (95% CI 7.5-9.63),respectively [19].A high HALP score (hemoglobin (g/L) × albumin (g/L) × lymphocytes (g/L)/platelets (g/L)) were significantly associated with lower OS (HR 0.707, 95% CI 0.503-0.995)and DFS (HR 0.671, 95% CI 0.491-0.916)[18].Unfortunately, these models are either based on bioinformatics analysis, which inevitably increases patient costs; Or, age-specific prognostic models fail to fully account for serologically specific factors and PD-L1 expression.Contrasted and these models, the advantages of NTRS are high advantage/cost and simple admittance to explicit markers in the determination factors.These assist clinicians with making a fast evaluation in light of the serological indicators of the patient.Apart from the above, the data in these studies [14-17, 21, 24-36] have only highlighted the value of tumor or nontumor-related biomarkers or models as single or independent factors for OS, PFS, therapeutic approaches and clinical responses in NSCLC patients, but have ignored the improvement in prognosis and prediction when combined with PD-L1 expression.Our study suggested that the risk score could more accurately predict PFS and OS in combination with PD-L1 expression in both the training and validation cohorts.This result has rapidly and easily complemented the application and reference significance of PD-L1 expression for clinical practice and has improved the predictive accuracy of OS and PFS.At present, PD-L1 is the main biomarker to judge the use of immunotherapy.Contrasted the above prognostic models, the combination of NTRS and PD-L1 expression can assess prognosis and response, instead of depending on a prognostic model alone.This is the main significance of the establishment of NTRS.
Importantly, several studies [37,38] demonstrated that predicting DCB was preferable in a number of scores/models that were built on strong economic conditions and methodologies.Among them, three genomic biomarkers, namely TMB, intratumoral heterogeneity and loss of heterozygosity in human leukocyte antigen, and the support vector machine (SVM)_poly method were used to construct a DCB prediction model.The DCB multi-feature model exhibited better predictive value with AUC = 0.77 (cohort1) and 0.78 (cohort2) [37].Another study [38] illustrated that a novel 18 mutated genes score 18 (TMS18) (HR 0.307, P < 0.001) had smaller HR and P value than PD-L1 expression (HR 0.403, P = 0.005) in survival analysis in 240 NSCLC patients.Moreover, the AUC of TMS18 was significantly higher than that of TMB, and the combination of TMS18 and PD-L1 improved the accuracy.The above studies were performed on the basis of TMB-related markers or genes, which made it more difficult to put into clinical practice.Our analysis suggested that the risk scores are more valuable than PD-L1 Furthermore, primary tumor LNM status [17], number and size of tumors [24], PD-L1 expression status [10,11], and ECOG [21,33,39] have been identified as important discriminatory variables in patients treated with ICIs.These studies have displayed the prognostic and predictive value of these factors, but few papers have examined the relative importance of "celebrity" markers in combination with other factors for survival.In the era of ICIs, the optimal cutoff value for assessing a prognostic factor may have changed.One possible explanation is related to the different disease burden of patients.A classification and regression tree model identified presence of LNMs (single or 2 vs. ≧ 3) as the most important prognostic factor among wild-KRAS patients; whereas, CEA levels (< 6 ng/mL, ≧ 6 ng/mL and < 360 ng/mL, and ≧ 360 ng/mL) were identified as the most important prognostic factor in mutation-KRAS patients [40].This study demonstrated that two separate factors or two levels of a biomarker can be used as the optimal cutoff value for prognostic assessment.However, these have been less studied in NSCLC.Of note, in patients stratified by risk score, primary tumor LNM constituted the only prognostic "node" for OS (Fig. 5); and SCC levels (" < 2.5 ng/mL" and "≧ 2.5 ng/mL") were the most important prognostic risk factors in patients stratified by PD-L1 expression (Fig. S7).The possible reason was that patients with N2 metastases have more extensive metastatic events; whereas, patients with N3 metastases were highlighted by LNMs and had relatively less distant tumor burdens.In contrast, patients with SCC ≧ 2.5 ng/mL and PD-L1 ≧ 50% had less PFS benefit.The result may be explained that these patients had more tumorigenic factors leading to worse prognosis, diluting the prognosis based on PD-L1 expression alone.
Interestingly, previous studies have observed that malignant pleural effusion [27] (HR 1.52, 95% CI 1.01-2.29,P = 0.043) were independent predictors of shorter PFS even though patients' PD-L1 TPS was high.Additionally, PD-L1 expression was highest in adrenals, livers and LNMs, and lower in bone and brain metastases [28].In our study, we observed higher PD-L1 expression in the LNMs and lower PD-L1 expression in the pleural effusions and adrenals, compared with higher risk scores in adrenals, livers and LNMs and lower in bones.These responses might be influenced by small sample sizes, multifocal metastases, and inconsistent  [41] demonstrated a 30% response rate to a single ICI in brain metastases in patients with driver-negative, high PD-L1 expression and low CNS disease burden.Our institution [26] has previously proved that high FARI was associated with poorer OS in patients with bone metastases (P = 0.03).Patients with brain metastases with PD-L1 expression > 1% had a longer OS, which was consistent with previous studies, but there was no prognostic difference between the score level and frequent metastatic sites.This may be due to the absence of correlation between the new risk score and metastatic sites, or the inability to effectively assess the metastatic characteristics of patients.The development of biomarkers or risk models remains an active area of investigation.
The three indexes in NTRS are ALP, CYFRA21-1 and HDL-C.Many studies have shown that elevated ALP is a robust adverse prognostic factor for survival in metastatic neuroendocrine tumors [42], spinal metastatic diseases [43], hepatic metastases [44], etc.These metastases predicted worse survival times.CYFRA21-1 has been associated with metastasis and poor prognosis [45].CYFRA21-1 levels had stronger association with metastases and survival status than CEA and NSE in all lung cancer patients [46].In this research, Cox regression analyses showed that only CYFRA21-1 was an independent prognostic factor in ADC and SCLC patient subtypes.Current literature suggests a positive association between cholesterol levels and the risk and progress of cancer development [47].This is principally because of the way that raised HDL-C gives the energy metabolism and membranous components expected by malignant growth cells to advance cancer proliferation and improvement.Consequently, the rationale is poor indicators that promotes tumor proliferation and metastasis behind selecting the specific factors.
Limitations of this study include the retrospective design, single-institution experiences, small samples, and restrictions of the electronic case system.Also, patients with ECOG ≧ 2 were treated with a single ICIs, similar to that of patients with PD-L1 expression ≧ 50%, which had a confounding prognostic and predictive implications.Second, existing studies on this topic are lacking, although the value of FARI as a prognostic biomarker in the assessment of patients with solid tumors receiving immunotherapy was explored.Third, diets were standardized prior to laboratory testing, but patients with combined insidious metabolic syndrome could not be completely excluded.Last, the limited time available led to incomplete sections of the survival data.
There are several strengths in this study.It's a new prognostic and predictive model and served as an independent poor prognostic factor for patients with advanced NSCLC undergoing ICIs based on PD-L1 expression or ECOG status as monotherapy or combination chemotherapy.And this

Conclusions
In summary, a simpler, faster and more accessible NTRS was developed compared to PD-L1 as a prognostic and predictive model for advanced NSCLC patients treated with ICIs compared to PD-L1.The score consisting of 3 non-tumorigenic blood factors may provide a universal approach to improve biomarkers in NSCLC patients applying only PD-L1 expression alone and will hopefully serve as useful feedback information for improvements for prediction and prognosis of ICIs in advanced NSCLC.
However, future prospective studies are required to confirm this conclusion due to the study's small samples.

Fig. 1
Fig. 1 Time-dependent ROC curve analysis to compare the effect of NTRS, PD-L1 expression, tumor size, ALP, CYFRA21-1, SCC, HDL-C, NLR, PLR, ALB, FIB, HCT, and LDH in predicting A, C overall survival and B, D progression-free survival in A, B training cohort and C, D validation cohort.The vertical axis shows the corresponding area under the ROC curve at different time points, and

Fig. 4
Fig. 4 Comparisons between PD-L1 expression and NTRS e in predicting response to ICIs in NSCLC in A, B, E, F, I, J training cohort and C, D, G, H, K, L validation cohort.Differential analyses of A, C NTRS and B, D PD-L1 expression among patients with durable clinical benefit.Differential analyses of E, G NTRS and F, H PD-L1 expression among PR/SD/PD patients.Number of patients with PR/

Fig. 5 A
Fig. 5 A Kaplan-Meier curve model comparing the hierarchical association of prediction factors associated with OS and PFS among patients with LNMs.OS overall survival, PFS progress-free survival, LNM lymph node metastases